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Abstract 

In this paper, we study those polynomials, orthogonal with respect to a particular 
weight, over the unioin of disjoint intervals, first introduced by N. I. Akhiezer, via a 
reformulation as a matrix factorization or Riemann-Hilbert problem. This approach 
complements the method proposed in a previous paper, that involves the construction 
of a certain meromorphic function on a hyperelliptic Riemann surface. The method de- 
scribed here is based on the general Riemann-Hilbert scheme of the theory of integrable 
systems and will enable us to derive, in a very strightforward way, the relevant system 
of Fuchsian differential equations for the polynomials and the associated system of the 
Schlesinger deformation equations for certain quantaties involving the corresponding 
recurrence coefficients. Both of these equations were obtained earlier by A. Magnus. In 
our approach, however, we are able to go beyond Magnus's results by actually solving 
the equations in terms of the Riemann 0-functions. We also show that the related 
Hankel determinant can be interpreted as the relevant r— function. 
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1 Introduction 

The Chebyshev polynomials are those monic polynomials characterised by the property that 
max|7r„(x)|, x G [—1, 1], is as small as possible. Indeed, it is also known that 7r„ is orthog- 
onal with respect to over [—1,1]. The polynomials vr^ — the Chebyshev polynomials 
of the first kind — which satisfy a constant coefficients three term recurrence relations, can 
be thought of as the "Hydrogen Atom" model of those polynomials orthogonal over [—1,1]. 
These play a fundamental role in the large n asymptotics of the Bernstein-Szego polynomi- 
als which are orthogonal with respect a "deformed" Chebyshev weight, p{x) — x^, over 
[— 1, 1], where p{x) is strictly positive, absolutely continuous and satisfies the Szego condition 



Many years ago N. 1. Akhiezer and also Yu. Ya. Tomchuk 12; considered a gen- 
eralization of the Chebyshev polynomials, where the interval of orthogonality is a union of 
disjoint intervals henceforth denoted as 



For comparison with those of Akhiezer, we assume here /?o = —1, and (3g+i = 1. For later 
convenience, when the end points become independent variables we shall adopt the conven- 
tion. 



EH 




E := iPo, ai) U aa) U ■ ■ ■ U {Pg, pg+i). 



(1.1) 



0^2, Po, Pi, Pg-^l 



{6i, 62, Sg+i, 6g+2, S2g+2)- 



(1.2) 



Let 




(1.3) 



be defined in the CP^ \ E. The multi-interval analog of the Chebyshev weight is 




(1.4) 
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and is obtained from the continuation w{z) to the top of the cut, E. The generahzed Cheby- 
shev or Akhiezer polynomials P„ are monic polynomials orthogonal with respect to 
i.e., 



where /i„ is the square of the norm. 

In the construction of the Bernstein-Szego asymptotics over E, for polynomials orthogonal 
with respect to the weight p{t)w^{t), where p is an absolutely continuous positive function, 
exact information on P„ would be required. This would entail the solution of the "Hydrogen 
Atom" problem in the multiple interval situation. In the case of two intervals, [—1, a] U [/?, 1], 
Pn was constructed by Akhiezer with an innovation which we would now recognise as the 
Baker-Akhiezer function, associated with the discrete Schrodinger equation, namely, the 
three term recurrence relations, where the degree of the polynomials n is the "coordinates" , 
and z is spectral variable. Akhiezer, based his construction on the conformal mapping of a 
doubly connected domain, with the aid of the Jacobian elliptic functions, as a demonstration 
for his students, the applications of elliptic functions j3]. It is not at all clear how the 
conformal mapping could be adepted to handle the situation when there are more then two 
intervals. In the early 1960 's, Akhiezer and also with Tomchuk published several very short 
and very deep papers regarding the Bernstein-Szego asymptotics. Akhiezer and Tomchuk 
gave a description of Pn and Qn (the second solution of the recurrence relations) with the 
aid of theory of Hyperelliptic integrals in terms of a cerian Abelian integral of the third kind. 
However, certain unknown points on Riemann surface appear in this respresentation, later 
circumvented in [H]. 

In a recent work of A. P. Magnus (6j, a general class of semi-classical orthogonal poly- 
nomials, which includes the Akhiezer polynomials Pn, was introduced and shown that these 
polynomials satisfy a certain system of linear Fuchsian equations. It was also demonstrated 
there that the recurrence coefficients, as functions of the natural parameters of the semi- 
classical weights, obey the nonlinear Schlesinger equations, i.e. the differential equations 
describing the isomonodromy deformations of the Fuchsian systems. 

In this paper we will study the Akhiezer polynomials P„ using the Riemann- Hilbert 
approach introduced in the theory of orthogonal polynomials in [7j. This will allow us to 
exploit the well-developed Riemann- Hilbert and algebro-geometric schemes of the Soliton 
theory [S], P, JO] - with certain important technical modifications though, and not only 
re-derive the previous results of j3] and JH] but also unite them in a single approach and 
produce further facts concerning the Akhiezer polynomials. Specifically, in addition to the 
derivation of Magnus's equations, we will solve them in terms of the multidimensional G- 
functions, and we will identify the corresponding Hankel determinant with the relevant r— 
function, i.e. with one of the central objects associated with an integrable system, in our 
case - with the Magnus- Schlesinger equation. It should also be mentioned that part of our 
- formulae, namely the ones describing the recurrence coefficients and the related Baker- 
Akhiezer function, reproduce the known expressions obtained in the late 70s (the works of 




(1.5) 
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I. Krichever, D. Mumford, S. Novikov, and M. Salle) for the finite-gap discrete Schrodinger 
operators which were then intensively studied in connection with the periodic Toda lattice 
(see the pioneering paper of H. Flaschka and D. McLaughlin JT] and also j8| and |^ for 
more on the history of the subject). 

We would like to think of our paper as a tribute to the pioneering works of N. I. Akhiezer 
which layed the foundation for the construction, in the 1970's of the algebro-geometric 
method in the theory of integrable systems, whose modern "Riemann-Hilbert" version we 
are using here. 

2 Riemann-Hilbert problem 

According to the classical theory of orthogonal polynomials the monic P^, (with Po = 1 and 
P_i = 0) and the polynomials of the second kind. 



of degree n — 1, are linearly independent solutions of the second order difference equation. 



Following the general scheme of [7] (see also JH]), let us introduced the 2x2 matrix 
Yn{z) be defined for n = 1,2,.. and z G C as follows: 




(2.1) 



ZVn{z) = Vn+l{z) + bn+lVn{z) + anVn-l{z). 



(2.2) 




(2.3) 



where 




—ITTWiZ 



(2.4) 



Proposition 1. The function Yn{z) satisfies the following conditions. 



RHl. Yn{z) is analytic in C\E 




RH3. Yn{z)z 



/, z 



oo. 
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where lAz^ denote a neighborhood of a point Zq. The matrix valued function Yn {z) is 
holomorphic in ^J z — [3j and bj is defined by the equation, 

w{z) = {z- P,r''\-{1 + 0{z - (3,)). (2.5) 

TT 

We shall also assume that the branch of ^J z — [3j is defined by the condition, 

< arg(z — Pj) < 2n, if j < g, and — n < aig^z — Pg+i) < vr. 
In addition, we assert that 

detfi^^)(/3,) = 1^0. (2.6) 

Proof. 

Using the basic properties of the Cauchy integrals and the Plemelj formulae we directly 
verify that Yn{z) satisfies RHl - RH2. To check property RH3 it is enough to note that 
because of the orthogonality condition il.^) . we have (cf. j3E]) 

r Pr^W^^^ = f:4^[ Pn{t)w,{t)t^dt 

Je z-t ^ ^^+1 Je 




To prove RH4 we observe that the matrix product, 

is bounded near (3j (the singular terms in the first column cancel out), and hence the function 
Yn^^\z) defined by equation RH4 is indeed holomorphic in ^ z — [3j. To complete the prove 
of the proposition we only need to establish equation ()2.6|) . To this end, we notice that we 
have already established RHl - RH4 but short of equation ()2.6|) . One can see, however, 
that RHl - RH4 already yield even stronger statement. Namely, we claim that 

deiYn{z) = l. (2.7) 

Indeed, the (scalar) function det Yn{z) is holomorphic in CP^ has no jumps across E and 
has removable singularities at the end points of E] moreover, it approaches 1 as z — > oo. 
By the Liouville theorem, equation ()2.7|) follows. Equation fl2.6|) is a direct consequence of 
equation (j2.7p . The proposition is proven. 

Remark 2.1 Equation (j2.7j) can be also derived by using the first line of (j2.3|) and the 
Christoeffel-Darbooux formula, 

(letYn[z) = / w+[t)dt 

"-n-l J E Z — t 
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n—1 ^ 



[ Kniz,t)w+it)dt = Po{z)ho = l, {2.t 
Je 



or from the recurrence relations, 

dety„(z) = -^(g„(2;)P„_i(2;) - P„(z)Q„_i(z)) = 1. (2.9) 

iln-l 

Remark 2.2 Yn{z) also depend on {5j : 1 < j < 2g + 2} . 

Proposition 2. Conditions RHl - RH4 defines the function Yn{z) uniquely. 
Proof. If Yn{z) is another function that satisfies RHl - RH4 then := Yn{z)Y^^{z) 

is holomorphic for z E CP^ \ : < j < (7 + 1}. Furthermore, for z E Up., 



Y-\z) = Y^'^^~\z), (2.10) 

where Yn^^'' ^(z) is holomorphic (see equation ()2.6|) !) in z — fij. This implies, 

X„(^) = 0(1), (2.11) 
which in turn implies X„(z) is holomorphic for z E CP^, and X„(z) = /, for all z E CP^. 

The conditions RHl - RH4 constitute the Riemann-Hilbert problem whose unique so- 
lution is given by equation (2.3), due to Proposition 1. 

The Riemann-Hilbert problem RHl - RH4 together with the equation 

P^{z) = (F„(^))n (2.12) 

will be used as an alternative definition of the Akhiezer polynomials. Notice also that the 
asymptotic condition RH3 can be extended to the full Laurent series, 

F„(.)= fj + f:^).-, \z\>l (2.13) 

and from (2.3) we have 

- I 

l/K-i -Pi{n) 



m,in) = ( fl^""^ . ) (2.14) 
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where pi{n) is the coefficient of z"- ^ of Pn{z). Taking into account the recurrence relations 
()2.2|) . we have, 

_ K 

and 

K.+i = Piin) -pi(n + 1), 



and the following relations supplementing ()2.12|1 

hn = (mi(n))^2 (2-15) 
a„ = (mi(n))^2 ("^i(^))2i (2-16) 
bn+i = (mi(n))^^ - (mi(n + 1))^^ . (2.17) 

Therefore, all the basic ingredients of the theory of polynomials Pn{z) (including the polyno- 
mials themselves) can be obtained directly from the solution Yn{z) of the Riemann-Hilbert 
problem. 

Remark 2.3 In the a prior setting of the Riemann-Hilbert problem RHl - RH4, the 
condition RH3 can be replaced by the following weaker one 



RH4. Y^{z) ( 

\ b, 







^z - 



0(1), zr^P„0<j<g + l. i2.U 



3 Differential Equations 

Having obtained equations ()2.12j) - ()2.17p which represent orthogonal polynomials Pn{z) and 
the corresponding norm and recurrence coefficients in terms of the solution Yn{z) of the 
Riemann-Hilbert problem RHl - RH4, we can now use the powerfull techniques of the 
Soliton theory. Specifically, in this and the two following sections we will apply a certain 
modification of the standard Zakharov-Shabat dressing method (see e.g. 1^) to obtain the 
relevant differential and difference equations for the Akhiezer polynomials. The modification 
needed is caused by the presence of the condition RH4. This condition indicates the relation 
of the problem under consideration to the theory of Fuchsian systems. Indeed, our derivations 
will be close to the Zakharov - Shabat scheme and to the constructions of the Jimbo-Miwa- 
Ueno monodromy theory (see also were both methods are unified in a single general 
Riemann-Hilbert formalism) . 

To describe the change of Yn{z) with respect to z for a fixed ra, it is advantageous to 
transform the Riemann-Hilbert problem satisfied by Yn{z) in to a form where jump matrix 
has constant entries. To this end, put 



$nW=>^n(.)U .l^i^A- (3-1) 
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A direct computation shows that 



<l>n,-{z) = <l'n,+ (z) 



1 -1 

-1 



(3.2) 



To specify the behavior of the new function near the end points of the set E let us observe 
that the new (constant !) jump matrix admits the following spectral representation, 



1 -1 
-1 



-1 
1 



where 



This implies that the function 



1 
2 -1 



v/J^ 

1 



satisfies the jump condition ()3.2j) in the neighborhood of Pj. Indeed, assuming z G Uj3. fl E, 
we find, 

' 0' 



[^f\z)]-^^^^'\z) = (v^) + 







1 



(v/7^_ 

1 



-1 
1 



1 -1 

-1 



Hence the matrix valued fucntion 



$„(z)[$(^^)(z)]-i 



has no jump accross E and therefore is holomorphic in the punctured neighborhood W/3 \{/3j}. 
Observe in addition that in the product. 



2m 
b, 











1 1 

2 2 

1 



1 



(3.3) 



the negative powers of z — fij cancel out. 

Therefore we conclude that the product ^n{z)[^^^^\z)]~^ is in fact holomorphic in the 
whole neighborhood Up.. Similar is also true for the matrix product 



$„(z)[$(°^-)(^)]-i^$„(^) 



1 



-1 



in the neighborhood Ua of the endpoint aj. Here we shall assume that the branch of ^Jz — aj 
is defined by the condition, 

—71 < a.Tg{z — aj) < 71. 
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In summary, $n(-2) solves the following Riemann-Hilbert problem: 
$1. is holomorphic for 2 G C \ 

-1 

n 







$4. M-) = $(f^)(.)(v^ .EW,, 

$5. $„(^) = <l^)(^)(vi^ _\), 



where <l'i^^'*(z) and ^l^^\z) are holomorphic in the neighborhoods of the points Pj and a 
respectively. Moreover, the matrices and ^l^^\aj) are invertible. In fact, 





J' 



and 



1 / TTt 

b, V 2 



a/ — 

Y 2 

— ni. n 

a, V 2 ^ 



where is defined by the equation (cf. 12. 5j) 

w{z) = {z- ajY/^aj^il + 0{z - aj)). (3.4) 

We want to emphasize, that unlike the case of the Y - Riemann-Hilbert problem, in the case 
of the $ - Riemann-Hilbert problem the left multipliers ^l!^^\z) and ^l^^\z) are holomorphic 
with respect to z. 

Remark 3.1 From $1 — $5 it follows (independent of 1)3.11) ) that 

det<l>„(^) = ^. (3.5) 
w[z) 



Consider now, the logarithmic derivative of $„(2) 



Aiz,n) := -^<l>-\z). (3.6) 

az 
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Since all the right matrix multipliers in the r.h.s. of $2 — $5 are constant matrices, A{z, n) 
enjoys the following properties: 

Al. A{z, n) is holomorphic for z e CP^ \ {aj, f3j}, 
'n 

A2. A{z,n)^— + Q ( ^2 ) ' ^ — ' °° 

A3. /l(^,„) = l4(A)(ft)-^^-|^i»)->(/J,)+0(l), z~P„ 

n'o 

I z — aj 



By virtue of the Liouville theorem, it follows that. 

Mi, 

z — (5j ' ^ z — a 



BJn) A A An) 



A{z,n) = J2^.+J2^. (3-7) 

j=0 j=l 1 



where 



^.W:=-^*^H«.)(J o)*^^-'K)--^^nK)(j o)nr'K)- (3.9) 



Note also. 



9+1 9 / n 

J]5,(n) + 5]A,(n)= 



Using (2.3) and RH4 give 



-n + 1 



^^^-^^l ^ b,P„_.(/.,)An- J ■ ^^-^^^ 

We conclude this section by recording the linear matrix differential equation with Fuchsian 
singularities at {Q;j,/9j}, mentioned in the Abstract, 

= A(^,n)$„(z), (3.11) 
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with A{z,n) defined by (3.7), (3.8) and (3.9). Furthermore, using the second hne of (2.3), 
the matrix valued residues are expressed in terms of the evalutions of the polynomials at the 
branch points: 

_ If QnmPn-im/K-i -QnmPnm \ ..... 

' 2 VQn-l(/5,)P„-l(/5,)//iLl -Qn-liP,)Pnif3j)/hn-lJ ' 

.. X _ 1/^ Pn{aj)Qn^l{aj)/hn^i -Qn{aj)Pn{aj) \ , . 

' 2\Qn^i{a,)Pn-M,)/hl_, -P„_i(a,)Q„(«,)//i„_ J • ^^''^^ 
Note that from (jSIHI) and (jSSl) it follows that 

trS,(n) = -^(Q„_i(/5,)P„(/?,) - Q„(/5,)P„_i(/5,)) = 1/2 
det Bj{n) = 0, 

and 

trAj(n) = -—^ — iQn{<yj)Pn^i{<yj) - Pn{oij)Qn-i{<yj)) = -1/2. 

detAj(n) = 0. 

We note that this leads to a discrete analogue of the "Wronskian" relation, 

Pn-l{z)Qniz) - Pn{z)Qn-l{z) = K-l, 

which, of course, can be independently derived from the recurrence relations. 

As it has already been mentioned in Introduction, equation (j3.1ip . even for more general 
weights of the type Y\_j{t — 5jY^, was first obtained in JU]. In 6J the Riemann-Hilbert problem 
is not used explicitely; rather, the author analyses directly the monodromy properties of the 
function Yn{z), i.e. the approach of [6] is based more on the ideas of Jl] than of [S]. It is also 
worth mentioning that our approach can be extended to the general semi-classical weights 
without any serious modifications. 

4 Derivatives with respect to the branch points. 

In this section we determine differentiation formulas for ^n{z) with respect to {a^, (3j}. First 
let us consider the logarithmic derivative of ^n{.z) with respect to a particular 13 j] 

VA^) := (4.1) 

and note that Vj{z) has the following properties 

VI. Vj{z) is holomorphic for z G C \ 
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V2. Vj{z) = 0{I/z), z-^oo. 





VS. v,iz) ~ + o(i). 

(4.2) 

By comparing with ()H.8|1 and again invoking the Liouville theorem, we conclude that 

VA^) = (4.3) 



which imphes 



A similar analysis gives 



z-pj 



d,M^) = -^^n{z). (4.4) 

Z Pj 



z — aj 



5 Difference Equation. 

Consider the "difference logarithmic derivative" 

f/„(z) := $„+i(z)$;i(^) = Yr,+^iz)Y~\z). 

Taking into account that all the right matrix multipliers in the r.h.s of RHl - RH4 are 
constant with respect to n we conclude that Un{z) is an entire function. Moreover, from the 
asymptotics (|2.13j) we have that 



£/„(,) = ( / + !!!l(!i±I)) -3 (; _ !!!;M) + o (i 



^(J o) + (o o)"(o ) + (- ) ' ^^oo- 



Appealing once again to the Liouville theorem, we conclude that Un{z) is linear function in 
z defined by the equations 



Uniz) 



z + (mi(n + l))ii - (mi(n))ii -(mi(n))i2 
(mi(n + l))2i 

i/K 
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where in the last equation we have taken into account ()2.15p - ()2.17|) . To summarize, the 
difference equation for the function ^n{z) reads 

Of course, equation ()5.H) is just the matrix form of the basic recurrence equation ()2.2|) . 
Nevertheless, we gave its "Riemann-Hilbert" derivation to emphasize the "master" role of 
the Riemann-Hilbert problem RHl - RH4 in our analysis. 



6 Schlesinger Equations and the Hankel Determinant. 

With the unified notation mentioned in the Introduction, we write, 

23+2 , . 

A{z,n) = J2^, (6.1) 

and the correspondence, 

(Ai(n), Ag{n), Bo{n), i?,+i(n)) ^ {C,{n), C,(n), C,+i(n), C2,+2(n)). (6.2) 
Note that, Cj{n), depend on 6jS. We of course have, 

23+2 / \ 

a,$„(^) = V^^$„(^), (6.3) 

dsM^) = -^^n{z). (6.4) 

Z — Ok 

Applying on (6.4) gives 

o o / N Ck(n) ^ Ck(n) '^•^ Cj(n) ^ 

S_,a^*4,) = _LJ_t„-_Aig_4J.j„, (6,5) 

and dsi^ on (6.3) gives, 

OsAM^) = T^^n + E '-^^n - ff (6.6) 

{z-6ky ^ z-6, yj^^ z-6j J z-6k 

Since d^ds^^^n = ds^^n and det $„ ^ 0, we get, 

E^^"^^ ds^Cjjn) _ [Cj(n),Cfc(?2)] 



^ (5j - 4 V^~'^i z-SkJ' 
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We now send z to a particular Sj in (6.7), with j ^ k, and find by equating residues, 

a,.C,(„) = ISM4M. (6.8) 

Oj — Ok 

If j = k, tlien a similar calculation gives, 

— Ol — Ok 

The equations (6.8) and (6.9) are the Schlesinger Equations satisfied by Cj{n). This is the 
equation first derived for the general semi-classical orthogonal polynomials in j6j. We are 
now going to move beyond the results of |6j and show that the corresponding r - function 
can be identified with the Hankel determinant associated with the weight w+(t). To this end 
we first recall Jimbo-Miwa-Ueno definition of the r -function. 
Let fi(^) be the one-form, 



dSj — d5k 
5i — 5k 

i<j<k<2g+2 



n'^'\5,,...,52g+2) := Yl tr(C,(n)Cfe(n)) 



J2 ti {Cj{n)Ck{n)) d\n\5j -5k\, (6.10) 

i<j<k<2g+2 

then it can be verified [H] using the Schlesinger Equations that, 

dfi(^) = 0, (6.11) 

which implies that, localy, fi'-^-* is an exact form. The r— function of the completely integrable 
system of partial differential equations (6.8) and (6.9) is then defined by the relation, 

=rflnr„(5i,..., 523+2). (6.12) 
In the theory orthogonal polynomials, the Hankel determinant, 

Dn[w+] := det ( ! t^^^w+{t)dt\ , (6.13) 

has two other equivalent expressions, 

Dn[w+] = ^ I ... I Y\_ {Xj - XkfYi^+^^l")^^^' 

l<3<k<n 1=1 

n-1 

= (6-14) 

j=0 
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It is to be expected from the structure of the Riemann-Hilbert formulation that, Dn, con- 
sidered as a function of is the r— function for this problem. To understand this, 
we require the derivatives of hn w.r.t. to 6^- To begin with, we use that, 

ds.Mz) = -^^<^n{z), (6.15) 



must be compatible with 



= ( oN^nW- (6-16) 



This results is 

z-bn+i -hn\ Ck{n) Ck{n+ 1) f z - bn+i -K 



hn J z — Sk z — Sk \ hn 

—dsi^bn+i —ds^hn 



-{1/ hn)dsi^ In hn 
which holds for all z e CP^ \ ...52g+2}- Putting 2; = oo in (6.17), gives, 

-Cf(n + 1) ; \-{l/hn)ds,\nhn 

which implies, amongst others, 

ds.hn = -Cf (n). (6.18) 

Lemma 1. Let the asymptotic expansion of A{z,n) about 2; = 00 be 

00 

A {z, n) = J2 A(n)^-'=-' (6.19) 

fe=0 

where 

2g+2 

Akin):=J2CAn)S^{n). (6.20) 

Then the first two Ak{n) are 

AW = {; i!„), (6.21) 



(6.17) 
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where Ci = ^1]=! iPj " "i)- 

Proof. Putting (3.1) into (3.11) we find 

Ay„(.-) + r„w(« _^;^^(^,)=A(.-,„)y„(.). (6.23) 

Expansion of (6.23) in z^^ gives the desired results. 

Theorem 2. The Hankel determinant is the r function of the Magnus - Schlesinger Equa- 
tions. 

Proof: We start by equating the residues of (6.17) at 2; = 6j. This gives, 



Ur,{Sj)Cj{n)^Cj{n + l)Ur,{Sj), 



or 



Cj(n + 1) = Ur,{5j)Cj{n)U-\5j), (6.24) 

where 



Un{z) 



l/hn 



A simple calculation shows that 



where a- :— ( ^ ^ ] ■ Now, 



U-\z)U^{z') = ^^'^ = I + '-^a., (6.26) 

1 0, 

dlnr^ = ^ ds^^ InTn dSk, (6.27) 

k 

where (from (6.10)), 

which leads to 

Tn+i _ V- tviCjin + l)C,in + 1) - C',(//)C',(n)) 



In 



'n+l _ 



r,, (5i — (^i- 
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^ triUniSj)Cjin)U-\5j )UMCk{n)U-\Sk) - C,(n)Cfc(n)) 



km - 

^ (5,- - Sk 

tr[(/ - ^a-)C,(n)(/ + ^a-)C.(n) - C,(n)C,(n)] 

'^i - '^"^ ■ i ■ ^ 



A calculation shows that the term [...] in (6.29) is 



We also note here some useful identities; 

tr(C,(n)a_Cfc(n) - a.Cj{n)Ck{n)) 

and 



= (n)(C^^(n) - (n)) - Cf{n){Cf{n) - Cf{n)), 



tr((7_C,(n)(7_Cfc(n)) = Cf (n)Cf (n). 

Therefore 

In ^ = 5^ {Cj'{n){Cl\n) - Cl\n)) - Cl\n){C]\n) - Cf{n))) 

- ^Y.^^-^k)Cf{n)Cl\n). (6.30) 



/)2 



To simplify the r.h.s. of (6.30) we note, from (6.20), (6.21) and (6.22) 

Y,{C]\n)-Cf{n)) = 2n-l, 

i 

Using these, and J2kmfk = -/j + Z^fe/fe, the r.h.s. of (6.30), becomes, 

Cfin) C]\n) C]\n) 

—y^ 2^(<^fe (^) - Cfe (^)) + 2^ 4Cfe (n) = — ^- — , 

" A; " k " 
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Finally, using (6.18) 



ds^ In ^ = ^ = ds^ In /i„. (6.31) 



Summing over n from to A^ — 1, we conclude that tn is a constant multiple of Dn, where the 
constant is independent of {6j}'^?^^. Since the r - function is defined up to such a constant, 
we can assume that the constant is unity; 

rjv(5i, ^29+2) = Dn[w+]. (6.32) 

Remark 6.1. It is worth mentioning that equations ()6.24|) follow also (by putting z = 6j) 
from the equation 

dll (z) 

Aiz, n + l)f/„(^) - U^iz)Aiz, n) = (6.33) 

which, in turn, is the compatibility condition of the basic Fuchsian equation ()3.11|) and the 
difference equation ()5.1|) . This is the matrix form of the so-called Freud equation which in 
principal can be written for any semi-classical polynomials - see ^| and (and also 7j). 
In the physical language this is the "discrete string equation" corresponding to the weight 
w+{t). More precisely, equation (j6.33p is the (discrete) Lax representation of the Freud 
equation which manifests its integrability from the algebraic point of view: linear equations 
fl3.11|) and ()5.H) form a Lax pair for the Freud equation (cf. jTj, [HI). 



7 Non-linear difference equations. 

As explained in Remark 6.1, the matrix equation ()6.33|) should lead to the nonlinear 
difference equations for the recurrence coefficients, following the genre of the Freud equations 
for the Akhiezer polynomials. To this end, we rewrite (6.24) elementwise, by first specializing 
6j to Uj and second to (3j. This will produce six difference equations, relating polynomial 
evaluations at the branch points and the recurrence coefficients. For later convenience we 
introduce four quantities 

ri"^ := XT^Pn(«,)Qn-l(«,), 
llln-l 

R^n^ ■■= ^^n(«,)Q„(«,), 

Thus by specializing to a^, Cj{n) becomes, 

R'nljK^i -rl"^-l/2 
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where we have taken onto account that the trace of the above is —1/2. In component form 
(6.24) is equivalent to, 

r^:l+r^:^ + \ = 4"H«.-&n+i) (7.1) 
a„+ii?Si - a„i?t\ = (6„+i - a,) (/?!")(&«+! - «,) + 2rW + 1) . (7.2) 

Note that out of the four possible equations, the 21 element is a tautology and the 11 and 
22 elements are equivalent. Similarly, specializing to j3j, Cj{n) becomes 



■ ^n-ll "-n-l ' n 

where the trace of the above is 1/2. In component form, (6.24) becomes, 

ri^l + ri^^ + \ = (/5,-&n+i) (7.3) 

a„+ii?l^i - a„i?i!.^ = (/5,. - 6„+i) (^i^H/?. - &n+i) - 2rf - . (7.4) 

In addition to these we have 

= ri")fr(") + i^ (7.5) 



2^ 

auR^^'R'^l, = fi+rfV (7.6) 



since det Cj(n) = 0. The equations (7.1) - (7.6) are the difference equations mentioned 
above. We should be able to eliminate, ri"\ rif\ and Rn'^ from these to obtain non- 
linear difference equations involving only a„ and These equations, are also discussed in 



8 The (Ji Riemann-Hilbert Problem. 

In this section we shall solve the Riemann-Hilbert problem RHl - RH4 for the Akhiezer 
polynomials in terms of the - functions. To this end we will need a further transformation 
of the Riemann-Hilbert problem satisfied by $n(^) to the so-called ai problem, first appeared 
in the theory of algebrogeometric solutions of integrable PDEs (see jTH|, PH]). 

We notice that since the matrices ^ ^ J ^ and cti have the same simple spectrum, 

they must be similar. Indeed, we have 
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Therefore, if we define 



° \ . . X A 1 



then the jump matrix of the new function becomes ai. The left diagonal constant matrix 
multiplier is introduced to normalize the asymptotic behavior of the function n{z) Sit z — oo: 

\ / n \ / \/27ri 



1 i I ) \ -Jf Ml -1 



i\\ / 





Taking also into account that 

we can reformulate the Riemann-Hilbert problem in terms of "ifni^), ^ follows. 

^n(-2) is holomorphic for 2; G C \ 

A oj /I -1 



= " 7 h (8-2) 



1 1 



1/2 

0^/^1-1 
1 1 



where is holomorphic in the neighbourhood of 2; = aj and det ^i°^^^(Q;j) 7^ 0,i.e., 



00 



fc=0 

20 



Similarly, '^n^\z) is holomorphic in the neighbourhood of z = [3j and det'^'^n^\[3j) ^ 0,i.e., 

oo 

= E ^nk\^ - det ^^^i^ ^ 0. 



k=0 

It is also worth noticing that the matrix products 



v/J^ 0\ (I -1 
and 

0\ A -1 



1/ 1 

have an exact ai— jump matrix in the respective neighborhoods. 

Remark 8.1. From \1/1 — \l/5 it follows (independent of ()8.1|l ) that 

det^„(z) = (8.3) 
7rti?(2;) 

Remark 8.2. The function \l/„(-2), in terms of Pn{z) and Qn(2;), is given as: 

/ ilTw{z)Pn{z) — Qn(z) i-Kw(z)Pn(z)-\-Qn{z) \ 

rr"! — — I ""(""^ 1 (R A'] 

\ h„-iw{z) h„^iw{z) / 

and all the properties listed in \1/1 — \l/5 can be deduced from this representation. It is worth 
emphasizing here that our approach does not require this formula. Our logic is: The initial 
Riemann-Hilbert Problem for Yn{z), quite generally posed, is transformed via (8.1) to the 
o"! problem which in turn leads to the equations (8.2) and (8.3) by the completely general 
principals of the Riemann-Hilbert problem. 

Let us now solve the cxi problem defined by \E'l — \l/5, however, without any reference 
to (8.4). The philosophy we adopt here is similar to that in the asymptotic analysis of 
orthogonal polynomials via the Riemann-Hilbert problem (cf. jT2], JB]): We simply "forget" 
the explicit formulas involving polynomials. 

Introduce the genus g Riemann surface 9^ defined by 

9 

y' = iz- Po){z - /3,+i) - «.)(^ - Pj), 

i=i 

and let \I/„(P), where P = {z,y) eDI be the vector Baker-Akhiezer function determined by 
the conditions: 

BAl. ^n(-P) is meromorphic on 9^ \ oo^ with the pole divisor, 

(C(P)) = -X^«, 
i=i 
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BA2. The behaviour of ^n(-P) at oo^ is specified by the equations, 

*„(/>) =((;)+o(i)).». p^oo+. 

in other words, 00+ is a pole of order n and oo~ is a zero of order n — 1. Here as usual, 00^ 
means 

P oo"^ <^==^ z ^ 00, y ^ ±z^^^. 
Let TT : 9^ ^ CP"*^ be the projection, 

n{P) = z, P={z,y), 
and * : !H — s> SH* be the involution, 

P^P* = (^,-y) if P={z,y). 
The main observation (cf. JH]; ^01) is that the matrix function, 

vl/„(z) := fvi;„(p),^„(p*)), (8.5) 



where vr(P) = 2;, and P 00^ as 2; — 00, solves the RH problem \E'l — \E'5. 

1. Indeed \E'l is satisfied by construction since ()8.5p defines "ifniz) uniquely as an analytic 
function on CP^ \ E. 

2. If z ^ E from the "+"-side (or from above the cut), then 

P^{z,y+{z)) = P^ 
P* ^iz,-y+{z)) = {z,y.{z))=P.. 

If z E from the "-" side, then 

P^{z,y4z)) = P_ 
P*^{z,-y4z)) = iz,y^iz))=P+. 



Hence, 



and it follows that, 



^^n+iz) = (^n(P+), ^n(P- 
^^n-iz) = ^n+(^)fTl, ZeE, 
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and therefore \l/2 is satisfied. 

3. We have by construction, z —>■ oo imphes P — * 00+ and P* — > oo~ . 
Therefore from BA2, 

*„w = (/ + o (!))(;" JL), (8.6) 

which shows that \I/3 is satisfied. 

4. The function is analytic in the neighborhood of P = /3j as a point of the Riemann 
surface IH. The local parameter at the point (3j is the square root of z — j3j. Therefore, in 
the neighborhood of P = (3j we have, 

00 

^n{P) = Y.^,,{z - P,f/\ (8.7) 

fc=0 

00 

^n{P*) = (8.8) 



k=0 



SO that 



k/2 



.fc=0 fc=0 



This in turn implies that the function \E'„ (z) defined by the equation \E'4 is a holomorphic 
function of z. Indeed we have 

¥^^^\z)^^4z)(^\ 1 

= 2 E [^^-^ - (-1) V.fc] - /5,)'''/', E I^:'-'^ + (-1) - /?,)'/' 

\fc=0 fc=0 

00 00 



J=0 1=0 

5. Since P = is a simple pole of ^^(P), the Taylor series (j8.7|) and (j8.8|) shoud be 



replaced by the Laurent series, 

00 

^„(P) = ^4(^-a,)'=/^ 

k=-l 

00 

#„(P*) = ^(-l)'=4(^-a,) 



A;=-l 



The rest of the arguments is literaly the same as in the /3 -case, and we have that the function 
'ifl^^\z) defined by the equation ^^5 is holomorphic ai z = aj. 



23 



Our final observation is that already established properties imply (8.3) (cf. our "Riemann- 
Hilbert" proof of ()2.7|) above) and hence the inequalities, 

det^^)(«,)^0, det^f)(/3,)^0. 

We now come to the 0— formula for \E'„(P). First we assemble here for this purpose 
some facts about the Riemann surface realized as a two-sheet covering of the z plane in 
the usual way and with the first homology basis depicted in the figure below. Let 



ai 




Figure 1: The dash curves represent the parts of the cannonical loops lying on the lower 
sheet. The lower (upper) sheet is fixed by the condition that it contains the point oo+ ( 
oo~). 



J 



be a set of normalised Abelian differentials of the first kind. As it is usual for a hyperelliptic 
curve, we shall chose the differentials dojj according to the equations, 

duj = 2^{A ^)jk dz, 

k=l y 

Ajk = / dz. 

Ja, y 

The invertability of the matrix A is a (relatively simple) classical result. We refer the reader 
to the monograph |Tn] for the basic general facts concerning the theory of functions on the 
Riemann surfaces (see also chapter 1 of POj)- Let us also introduce the normalized Abelian 
differential of the third kind, having its only poles at 00=*", 

MP) = " + + - + 

y 

with vanishing a— period; 

dn = 0, j = l,...,g. 
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The above g conditions uniquely determine [12] the coefficients, {Xj}j^l. Put 

n{P) = I dn. 

One easily deduces, 



fi(P) = ± ( Inz-lnC(E) + ( - ) ) , P^oo^, (8.9) 



where 



C{E) = exp I - / I ^/=° ' --\dz\. (8.10) 



y(2;) z 



(We recall that = 1.) Finally, the Riemann B— function of g— complex variables s & C^, 
is defined with the aid of the period matrix 

Bjk '■= / dujj, 

as follows: 

e(s) = e(s; B) ■= ^ exp ( m{t, Bt) + 2-Ki{t, s) ) . 

Here are the fundamental periodic property of the G— function: 

e(s + n + Bm) = e— ^(^^.'^)-2-iW^)e(^, (8.11) 

and the obvious symmetry relation: 

ei-s) = e{s). 



Observe now that BAl — BA2 imply the following properties on the components of 

^n(P). 

\l/„i(P) is meromorphic on 91 \ {oo^, oo^} 

9 

^„i(P) =z" + 0(^"-^), P^oo+ 

^„i(P) =0(^-"), P^oo-. (8.12) 
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Similary for '^n2{P), 

vl/„2(P) = 0(z"-i), P^oo+. (8.13) 
By standard technique of the algebrogeometric method ( see e.g. [ID]), we get, 



^nl(P) = 
^n2(P) = 

where 




and Cj form the vector of the Riemann constants (see again JH] and [ID])- Indeed, by the 
Riemann theorem (see e.g. [19]), the first 0— functions in the denominators has zeros exactly 
at the points aj; the front exponential factors provide the needed asymptotic behavior at 
oo^; the first 6— functions in the numerators, by virture of the periodicity property ()8.11|) . 
ensure the single-valuedness; and, finally, the P-independent 0-factors together with the 
back exponential factors provide the needed normalizations ad oo^ (cf. fl8.12p and ()8.13|) ). 
We also assume that we choose the same path between and P for all the integrals 
involved^ . 

The formulae above can be simplified. To this end we observe that 

where the path of integration from to lies on the upper plane of the upper sheet. 
Therefore, moduli the lattice periods, 

9 

D, = 1 + ^P^.,((^-A: + 1). 

fc=i 

^Alternatively, one can choose for each integral its own path. In this case though the paths must not 
intersect the basic cycles. 
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In other words, the vector D belongs to the latice Z^ + i?Z^ and hence (property (|8.11|) again) 
can be droped from the above formulae for \E'„(P). This yields the following simplified 0— 
representation for \E'„(P). 

^ni(P) = e-^(^)^%^- ^ , 7-"^ \. C-{E), (8.15) 

ef/J rfcu + (n-i)L) do) 

e{Q^du-in-l)L) 

We conclude the G— function solution of the Akhiezer Riemann-Hilbert problem by 
noticing the following equation for the vector L of the b - periods of the integral Q{P). 

+ POO" POO'^ 

du+ dQ = -2 I du, (8.16) 

The equation is just the classical Riemann bilinear identity (see e.g. JP] or ^^) applied to 
the pair of the Abelian integrals cJ(P) and VL{P). 

Remark 8.3 Using equation ()8.14|) . one can check directly, with the help of the periodic 
condition ()8.11|) . that the theta function. 




has the points aj as its zeros. 

Remark 8.4 The reader should not be misled by the formal possibility to diagonalize 
simultaneously the jump matrices of the Riemann-Hilbert problem ^^1 — ^^5 (which all are 
equal to ui) and by apparently following from this conclusion that the problem can be 
reduced to the scalar one on the complex plane and hence solved without any use of the 
O— functions. The obstractions come from the end points a^-, (3j and from the point at 
infinity, where the function \E'„(^) must have the singularities specified by equations ^5, ^4 
and \E'3, respectively. These singularities can be alternatively discribed as the addition jump 
conditions posed on the small circles around the end points and on the big circle around the 
infinity. The relevant jump matrices are 

l/yj^ 0\ A -l\ f^^^, 0\ A -1 
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and 



z 



n 









n + 1 



) 



respectively. Posed in this form, the o"i Riemann-Hilbert problem becomes the regular one 

- no singularities different from the jumps are prescribed. At the same, the additional jump 
matrices depend on z and the whole new set of jump matrices can not be simultaneously 
diagonalized. The only way to circumvent this obstacles, and not to use the 0— functions, 
is the equation (8.4) which indeed gives an explicit representation of the solution of the cxi 
Riemann-Hilbert problem in terms of the elementary functions and their contour integrals. 
The 0— function representation ()8.15|) for the solution \l/„(2) obtained in this chapter has an 
important advantage comparing to (8.4). It reveals the nature of the dependence of '^n{,z), 
and hence of the Akhiezer polynomials themselves (see ()9.1|) below), on the number n, as n 
varies over the whole range 1 < n < oo (see jHj for more on the use of the B - representations 
in the analysis of the Akhiezer polynomials). Simultaneously, the comparison of equations 
(8.4) and ()8.15|) might, perhaps, be used to derive some new nontrivial identities for the 
hyperelliptic 6— functions. 

Remark 8.5 Up to a trivial diagonal gauge transformation, the matrix function '^n{z) 
satisfies the same Fuchsian equation that is satisfied by the function <^n{z). Note 

that the corresponding monodromy group is very simple; indeed, it has just one generator 

- the matrix o"i. Once again, the reader might be wondering about the appearance of the 
highly nontrivial theta-functional formulae in the describtion of the function \l/„(z) which 
gives the solution of the corresponding inverse monodromy problem. Similar to the previous 
remark, the explanation comes from the fact that the solution \l/„(2;), in addition to the 
given monodromy group, must exhibit the local behavior at the singular points indicated by 
the conditions \l/3 — \l/5. This situation is typical in the theory of the finite-gap solutions of 
integrable PDEs ^ (see e.g. [Hj and [TU]). 

9 A list of the G - formulae. 

In this section, we give formulae expressing the polynomial Pn{z)-, recurrence coefficients a„, 
hni the square of the weighted norm /i„ and the Hankel determinant in terms of the 6— 
functions. The expressions will be derived as simple corollaries of the equations ()8.5|) and 
(|8.15p representing the solution \l/„(-z) of the Riemann-Hilbert problem — \E'5 in terms of 
the - functions. 

^Another example of an apparently simple but nontrivialy solved invesre monodromy problem can be also 
found in the theory of integrable PDEs. It is provided by the multi-soliton Baker- Akhiezer function whose 
monodromy group is just trivial. Of course, the formulae in this case are simplier than the finite-gap ones - 
they do not contain the O— functions. At the same time, the answer is still rather complicated; in fact, it 
involves degenerated functions corresponding to the singular curves of genus zero. 
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From ()8.1|) it follows that (see also 

P^{z) = (F„(^))n = (^„(^))ll + {^n{z))l2. 

This together with ()8.5p and ()8.15|) leads to the following B - representation of the Akhiezer 
polynomials, 



PJz) 



Q {^^ + du) e-"(-) + e{nL- j;^^^ duj) e- 



■nQ(z) 



X . (9.1) 



where all the hyperelliptic integrals are taken in the upper sheet of the curve 91 (and along 
the same path). 

Remark 9.1 It is a simple but an instructive exercise to check directly, using equation 
()8.14|) . the similar equation for the integral ^{P), i.e. 

k 

and, once again, the periodicity property of the 6-function, that the right side of (j9.ip is 
indeed a polynomial. 

To evaluate the quantities a„, 6„, and hn we shall use the relation 



between the first matrix coefficients, ipi and mi, of the Laurent series 

^n(^)= / + $^^ , , k|>l, 



fc=l 



and 

oo 



respectively. In ()9.2p . the parameter k is defined via the expansion. 



w{z) = — 1 + - 

TTZ \ Z 
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and the matrix, 

' 



-1 oy ' 

should be added to the r.h.s if n = 1. Combaining equation ()9.2|) with the formula ()2.15|) we 
obtain that 

hn = 2(V'i(n))i2- 

On the other hand, let us introduce the coefficient matrix cjk, j,k = 1,2 hj the relations (cf. 
(pmi) and dHm), 

^„^(P) =;2« + cnz'^-i + 0(z"-2), P^oo+ (9.3) 
=ci2^-" + 0(z-"-i), P^oo", (9.4) 

and 

*n2(P) = ^-"+'+C222-" + 0(2-""l), O^OO-, (9.5) 
*n2(P) = C2iZ"-^ + 0(z"-2), P^OO+. (9.6) 

Then, it is obvious that 

iMn))jk = Cjk, (9.7) 
and, in particular, we arrive to the equation 

hn = 2ci2. (9.8) 

The coefficient cu, in its turn, can be immediately evaluated from the B— formula ()8.15p by 
letting P —>■ oo~ . In fact, we have 



9 ( duj — riL 
cu = C'-{E)^^ ^. (9.9) 



Taking into account the Riemann bilinear relation ()8.16p we can present the formula for hn 
in the following final form. 



e((n + i) l) 

hn = 2C2n(^) V n=l,2,..., (9.10) 



Q[{n-\)L 
ho := 1. 

An important direct consequence of this equation is the explicit G— functional representation 
for determinant of the (n + 1) x {n + 1) Hankel matrix: 

e((n + l)L) 

Dn+i[w+] = llh, = 2"(C(E))"("+i)^ ^ 



j=0 
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e(iL 



2"(C(E))"("+^) 



©((2- 



(9.11) 



A similar use of the remaining equations in ()9.2p . ()9.7|1 and the formulae ()2.1()j) . ()2.17j] 

leads at once to the 0— representations of the recurrence coefficients a„ and 6„: 



if n 



(9.12) 



e2((„-i)L) 



if n > 1 



and 



1 ^ 



e^((n^i)£) e;.((n-3)i;) ^ e^(lL) 

e((n-i)£) e((n~3)L) e(^L) 



(9.13) 



Here, 



0' := 



99 (s) 



Equations dnH), (j^TTIUl . dHUIl), (f^:!^ and were previously obtained in |3] by a direct 

analysis of Akhiezer's function defined as the sum (^q^ (j8.4j) ). In [H] it was 

also shown that the above formulae allow to identify the quantity C{E) as the transfinite 
diameter of the set E. We remind that in our approach, C{E) appears as a first nontrivial 
coefficient in the asymptotic expansion of the Abelian integral fi(P), see fl8.9|l and ()8.1()|1 . 
Finally, we should note that equations ()8.15|) . ()9.12|) and ()9.13|) . as the equations describing 
the eigenf unctions and the coefficients of a finite-gap discrete Schrdinger operator, have 
already been known ( see e.g. [20j) in the theory of the periodic Toda lattice. 
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